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^ : ABSTRACT 

^. 

^_\ ■ We present the newly improved Batalin-Fradkin-Tyutin (BFT) Hamiltonian formalism and 

CN ■ the generalization to the Lagrangian formulation, which provide the much more simple and 

f — I transparent insight to the usual BFT method, with application to the non-Abelian Proca model 

- ' which has been an difficult problem in the usual BFT method. The infinite terms of the 

effectively first class constraints can be made to be the regular power series forms by ingenious 

^. choice of Xa/3 and cj^^-matrices. In this new method, the first class Hamiltonian, which also 

,-C3 . needs infinite correction terms is obtained simply by replacing the original variables in the 

original Hamiltonian with the BFT physical variables. Remarkably all the infinite correction 

rS I terms can be expressed in the compact exponential form. We also show that in our model 

5^ I the Poisson brackets of the BFT physical variables in the extended phase space are the same 

structure as the Dirac brackets of the original phase space variables. With the help of both our 

newly developed Lagrangian formulation and Hamilton's equations of motion, we obtain the 

desired classical Lagrangian corresponding to the first class Hamiltonian which can be reduced 

to the generalized Stiickelberg Lagrangian which is non-trivial conjecture in our infinitely many 

terms involved in Hamiltonian and Lagrangian. 
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I. Introduction 

The Dirac method has been widely used in the Hamiltonian formahsm [1] to quantize the 
first and the second class constraint systems generally, which do and do not form a closed con- 
straint algebra in Poisson brackets, respectively. However, since the resulting Dirac brackets 
are generally field-dependent and nonlocal, and have a serious ordering problem, the quanti- 
zation is under unfavorable circumstances because of, essentially, the difficulty in finding the 
canonically conjugate pairs. On the other hand, the quantization of first class constraint sys- 
tems established by Batalin, Fradkin, and Vilkovisky (BFV) [2, 3], which does not have the 
previously noted problems of the Dirac method from the start, has been well appreciated in 
a gauge invariant manner with preserving Becci-Rouet-Stora-Tyutin (BRST) symmetry [4, 5]. 
After their works, this procedure has been generalized to include the second class constraints 
by Batalin, Fradkin, and Tyutin (BFT) [6, 7] in the canonical formalism, and applied to various 
models ^^^^ obtaining the Wess-Zumino (WZ) actions [11, 12]. 

Recently, the BFT Hamiltonian method [7] has been systematically applied to several linear 
[13] or non-linear [14] second class constraint systems producing the interesting result to find 
the new type of the WZ action which cannot be obtained in the usual path-integral framework. 
It is interesting to note that even the non-linear system requiring infinite iterations can be 
exactly solved in this method by some ingenious choice of the arbitrariness in the definition 
of {Xai3,uj°'^). All these works are based on the systematic construction of the first class 
Hamiltonian as a solution of the strongly involutive relation with the effectively first constraints. 
However, this method may not be easy one to much more complicated systems due to the 
complexity of the routine procedure [15]. 

In this respect we have recently suggested the improved method finding the first class 
Hamiltonian much simply and more transparently than the usual method and applied to several 
simple models, i.e., Abelian Proca model [16] and Abelian chiral Schwinger model [17]. However, 
these models do not show dramatically the power of our improved method due to the simplicity 
of the system. In this sense it would be worth considering the non-Abelian Proca model [18] 
which has not been completely solved in the usual method [8, 15]. This is one issue which will be 
attacked in this paper. However, even in this improved method the corresponding Lagrangian 
formulation was unclear due to the appearances of the time derivatives in the Lagrangian [16, 
17]. This is another issue on the formalism which will be also tackled in this paper. 

In the present paper, we present the newly improved BFT Hamiltonian formalism and 
also the generalization to the Lagrangian formulation, which provide the much more simple 
and transparent insight to the usual BFT method, with application to the non-Abelian Proca 
model. The model is most non-trivial model due to the necessity of infinite correction terms 
involving newly introduced auxiliary fields to perform the BFT's conversion of the weakly 



second class constraints into the strongly first class constraints. However this infinite terms 
can be made to be regular power series forms by ingenious choice of the X^p and cu^^-matrices 
such that our formulation can be also applied in this case. 

In Sec. II, we apply the usual BFT formalism [7] to the non-Abelian Proca model in order to 
convert the weakly second class constraint system into a strongly first class one by introducing 
new auxiliary fields. We find that the effectively first class constraints needs the infinite terms 
involving the auxiliary fields which being regular power series forms in our ingenious choice. 

In Sec. Ill, according to our new method, the first class Hamiltonian, which also needs 
infinite correction terms, is obtained simply by replacing the original variables in the original 
Hamiltonian with the BFT physical variables defined in the extended phase space. The effec- 
tively first class constraints can be also understood similarly. It is proved that the infinite terms 
in the BFT physical variables and hence first class constraints and the first class Hamiltonian 
can be expressed by the compact exponential form. We also show that in our model the Poisson 
brackets of the BFT physical variables fields in the extended phase space are the same as the 
Dirac brackets of the phase space variables in the original second class constraint system only 
with replacing the original phase variables by the BFT physical variables. 

In Sec. IV, we develop the Lagrangian formulation to obtain the first class constraints which 
complement and provide much more transparent insight to the Hamiltonian formulation. Based 
on this Lagrangian formulation we directly obtain the Stiickelberg Lagrangian as the classical 
Lagrangian corresponding to the first class Hamiltonian which is most non-trivial conjecture 
related to Dirac's conjecture [1] in our infinitely many terms involved model. This is also 
confirmed by considering the Hamilton's equations of motion with some modification term in 
Hamiltonian which is proportional to the first class constraints Gf . This disproves the recent 
argument on the in-equivalence of the Stiickelberg and BFT formalism. 

Sec. V is devoted to our conclusion and several comments 

II. Conversion from second to first class constraints 

Now, we first apply the usual BFT formalism which assumes the Abelian conversion of the 
second class constraint of the original system [7] to the non-Abelian Proca model of the massive 
photon in four dimensions [8, 15, 18, 19], whose dynamics are described by 



_ pa pa,nu _| , 

4^^ 2' 



S= I d'x [--F;,F'^'^^ + -m'A^^A^''^], (1) 



where F^, = d^Al - d^A^ + gT^A'^Al, and g^, = diag{+, -, 
The canonical momenta of gauge fields are given by 

6S 
< = —. — ^ 0, 



< = -^ = K (2) 

with the Poisson algebra {A"'^(a;), 7r|^(|/)} = (5^5ab5'^(x — y). The weak equahty ' ~ ' means 

that the equahty is not apphed before aU involved calculations are finished. [1] In contrast, the 

strong equality ' = ' means that the equality can be applied at all the steps of the calculations. 

Then, O" = ttq ^ is a primary constraint [1]. On the other hand, the total Hamiltonian is 



Ht = H^+ d^xu^'Ql (3) 

with the Lagrangian multipliers -u" and the canonical Hamiltonian 



Hr = / d^x 



-«)2 + -F^^F"''^ + -m2{(A"'°)2 + (A"'^)2} - AlQl 



(4) 



where 02 is the Gauss' law constraints, which come from the time evolution of 0" with Ht 
defined by 



= (DV,)" + m^A"'^ ^ 0. (5) 

Note that the time evolution of the Gauss' law constraints with Ht generates no more additional 
constraints if we choose only determines the multipliers as u"" ~ — (9j/l"'*. As a result, the full 
constraints of this model are 0q, {a = 1, 2) which satisfy the second class constraint algebra as 
follows 

A^^,(x, y) ^ {0^(x), 0j(y)} = (^ ^2% 'Jllc^c ) 5'(x - y), (6) 

where we denote x = (t, x) and three-space vector x = {x^,x'^,x^) and ei2 = — e2i = 1, and 
Q'^ = {D'^TTiY. Here we note that it is dangerous to use the constraints 02 ~ of Eq. (5) at 
this stage like as in Ref. [15] since in our case 0^ ^ of Eq. (9) only can be imposed. 

Now, we introduce new auxiliary fields $" to convert the second class constraints 0^ into 
first class ones in the extended phase space with the fundamental Poisson algebra 

{^,$^'^} = 0, 
{C(a:),<(y)} = ^.t(a:,l/) = -u;i"{y,x), (7) 

where JF = (A^, tt^). Here, the constancy, i.e., the field independence of uj"l^{x, y), is considered 
for simplicity. 

According to the usual BFT method [7], the modified constraints 0^ with the property 

{0^,0j} = o, (8) 



which is called the Abelian conversion, which is rank-0, of the second class constraint (6) are 
generally given by 

oo 

0^[^;$] = 0^ + ^e^(") = o; e(r)~($r (9) 

n=l 

satisfying the boundary conditions, O^[jF;0] = O^. Note that the modified constraints O^ 
become strongly zero by introducing the auxiliary fields <l>^, i.e., enlarging the phase space, 
while the original constraints O^ are weakly zero. As will be shown later, essentially due to 
this property, the result of the Dirac formalism can be easily read off from the BFT formalism. 
The first order correction terms in the infinite series [7] are simply given by 

Qf\x) = Jd'yX:1,{x,y)<^^,{y), (10) 

and the first class constraint algebra (8) of 0^ requires the following relation 

A^^^ix, y) + j d'w dh X- (x, w)u;2!{w, z)X%{y, z) = 0. (11) 

However, as was emphasized in Refs. [13 - 15], there is a natural arbitrariness in choosing the 
matrices cu"^ and X^^ from Eqs. (7) and (10), which corresponds to canonical transformation 
in the extended phase space [6, 7]. Here we note that Eq. (11) can not be considered as the 
matrix multiplication exactly unless X^j{y, z) has some symmetry, X^'^{y, z) = kX^p{z, y) with 
constant k because of the form of the last two product of the matrices / d^ztoj^ {to, z)XJ^'^{y, z) in 
the right hand side of Eq. (11). Thus, using this arbitrariness we can take the simple solutions 
without any loss of generality, which are compatible with Eqs. (7) and (11) as 



^:S{^.y) - 


= e^^6^'6'{^-y). 


K'p{^,y) - 


( -6"-^ 



'^'(x-y), (12) 

i.e., antisymmetric uj^^{x,y) and Xij{x,y) such that Eq. (11) is the form of the matrix multi- 
plication exactly [10, 13, 15] 

A':^Jx,y) - f d'u;dhXZ{x,u;)u;fAu;,z)XfJz,y) = 0. 



Note that X'^Jx,y) needs not be generally antisymmetric, while uj"l^{x,y) is always antisym- 
metric by definition of Eq. (7). However, the symmetricity or antisymmetricity of X^^(a:, y) 
is, by experience, a powerful property for the solvability of (9) with finite iteration [13] or with 
infinite regular iterations [14]. Actually this solution provides the latter case in our model as 
will be shown later. 



Now with this proper choice, the first order corrections of the modified constraints become 



ef ) = 


= -m^l 


ef ) = 


= m<!>l + £-r'^<!>ln^ 


= 


2m 



(13) 

where the cross operation 'x' represents {A x 5)" = f"-'^'^A^B^. Note that ©2 have the 
non-Abehan correction terms contrast to 0" . 

The higher order iteration terms, by omitting the spatial coordinates for simphcity, [7] 

1 

n + 2 



Qt^'^ = -^rn^^'A;x2^B',t^ {n > 1) (14) 



with 

n 71—2 



m=0 171=0 



are found to be non- vanishing contrast to the Abehan model [15, 16]. Here, cu^^ and Xjj: are 

dA dB dA dB 



the inverse of uj^J. and Ar;^«, and the Poisson brackets including the subscripts are defined by 



{^)^}(Q,P) 

r FiA f)R r)4 r)R 1 

(16) 



a,l3 L 



dQ dP dP dQ ' 

dA dB dA dB 



(9$° d<l>(^ (9$/5 5$o 



where {Q,P) and ($",$'^) are the conjugate pairs of the original and auxiliary fields, respec- 
tively. After some calculations it is not difficult to prove that 

Qa(n+1) ^ Q^ 

1 / ^ \ "+1 



using the mathematical induction. Here 'n-fold' means that the number of the parenthesis 
and the bracket is n. Note that all the higher order corrections 02 are essentially the 

non-Abelian effect. 

Hence in our model with the proper choice of Eq. (12) the effectively first class constraints 
to all orders become 



.2 

oo 



0^ = 0?-m<l>^, 

0^ = e^ + m$i + E 7-^(^1 [$^x($^x(...($^xfi)...))]J,.foH)- (1^ 



T'l {n + 1)! Vm 



with [$2 X ($2 X (■ ■ • ($2 X fi) ■ ■ ■))]o-fold - ^"• 

III. Physical variables, first class Hamiltonian, and Dirac brackets 

Now, corresponding to the original variables JF, the physical variables, called BFT physical 
variables in the extended phase space, within the Abelian conversion, JF, which are strongly 
involutive, i.e., 

{e^,^} = o (19) 

can be generally found as 

^[jr. $] = jr + ^ ^(«)^ ^W _ ($)« (20) 

n=l 

satisfying the boundary conditions, JF"''^ [JF; 0] = JF. Here, the first order iteration terms which 
are given by the formula 



^(i) = _$%^xT"^{e5,^} 



(^) 



f211 



become as follows 



If) = -$i-/^($^xfi)«. 



~a(l) 9 /^2 ^ ^\a 

m 

~a(l) _ Qa(l)_ ^22) 

The remaining higher order iteration terms which are given by general formula [7] 

^(n+i) ^ ^$/3^ x^^(G^)i") (23) 

n + 1 

with 

n n— 2 

l'-^FJ5 - Z^iy^ >-^ I(^)+Z^|y5 5->^ l('J') + |y5 .-^^''JC*), (24) 

m=0 m=0 



are also found to be as follows after some routine calculations 

n+1 I ( g'^ "+^ 



A 



a{n+l) 




{n + 2)\m^\mj L V V V ; ^^J(n+l)-fold' 



a: 



a(n+l) 



n + 



9V^\^2 



IJIUJ [$^x($^x(...($^xA.)...))]J^^^)_f^l, 



1 1 [9 



[n 



n+l 



nioUJ l^'x(^'x(---(^'x^^^')---))]n-fold' 



TT, 



a{n+l) 



TT, 



a{n+l) 




-fold' 



0^"+') = 0. (25) 

Hence, the BFT physical variables in the extended phase space are finally found to be 

1 fg 



^0 



A' 



^s+-*;-E 



n 



\{n + 1)! m? \m 



[$'x ($'x (•■■($'xfi) •■■))];_ 



fold' 



°° ^ f9 ^ r^2„.^2 



n=l 



n! vm 



^^ + ErT z: [$^x($2x(...($2xA.)...))]^. 



fold 



- 1 1^^ (£)" [^^ X ($^ X (. . . ($^ X 5.$^) . . ■))](,.i).fold' 



n=l 



1 

n! \m 



-EA - [^'x($ 



X 



(<^'x9,$2)...))]Lr 



-fold 



+ 



TT. 






^ 



TTn 



expj — ^^XJTTj 

m 



< + E-tP- [$^x($^x(...($^xvr.)...))];. 



fold 



(26) 



by defining the exponential of the x -operation as follows 



exp{Ax}B = B + AxB + -{Ax{AxB)) 



B+J2-[Ax{Ax{...{AxB)-- ■))]n-fold- 



n=l 



n! 



(27) 



Here, we note that the appearance of the infinite terms of correction are independent on the 
(non-Abelian) gauge group such that a genuine property of the Abelian conversion of the non- 
Abelian model of Proca theory. However in this form the gauge transformation property of 
the BFT physical variables is not so transparent. Hence we needs to consider more compact 



8 



expressions to achieve this goal. To this end, we, first of all, note that the following relation 

[i0, [i0, [■ ■ ■ [#, M] ■ ■ -llln-fold 

= (-1)"[0 X (0 X (■ ■ ■ (0 X M) ■ ■ ■))]n-fold (28) 

is satisfied such that the following useful formula are produced as 

CO 1 

e^^Me-^^ = E3[^</'J#J---[^</'>^]-- -En-fold 

oo (_l\n 

= E^['^x(</>x(---(</>x^)---))]n-fold 

n=0 "'■ 

= exp{-0x}M, (29) 



da ;^o(^ + l)' da 



J n- fold 



oo (_l\n Q 

- E T^\^ X (0 X (. . . (0 X ^-0) . . .))]n.fold' (30) 



n=0 

^bl •„ fabcrc 



where = (j)'^T°',M = M'^T^, [T^-^T ] = igf T'^ for any non-Abelian gauge group Q with 
generator T". Then, by using these formulas (29) and (30) A" and vf" can be re-expressed more 
compactly as follows 

A^[T;^] = --(diW-W-'y + {WA,W-')\ (31) 

ntiJ';^] = {Wtt^W-'Y (32) 

with the matrix valued factor W = exp{ ——^'^}. Furthermore, it is not difficult to show the 
following formula, by using the form (31) 

= Ft^in (33) 

Actually, these compact expressions of (31)-(33) can be easily expected ones in the BFT for- 
malism if we remind the spirit of the formalism itself. In this formulation, the BFT physical 
variables JF satisfy the Eq. (19), i.e., 

{e^,^} = 0. (34) 

But, if we remember that the first class constraints (here 0J^ contrast to second class constraints 
O^) can be generators of gauge transformations, it is clear that JF should be gauge invariant 
one in this BFT formalism. Actually, under the usual non-Abelian gauge transformation 

A^ - < = -'- {d,u . u-'Y + {uA.u-^y , 

F;.-^F;J={uF,,U-'Y (35) 



the combinations of (31)-(32) are the only possible combinations of manifest gauge invariance 
with the factor W transformation 

W ^W = WU-\ (36) 

This is the story of the BFT physical variables JF corresponding to the original phase space 
variables J-'. On the other hand, it is interesting to note that the BFT physical variables for 
the auxiliary fields, which were not existed originally, become identically zero as it should be, 
i.e., 

K = 0. (37) 

Similar to these BFT physical variables corresponding to the fundamental variables JF, the 
BFT physical quantities for more complicated physically interesting quantities can be found by 
the systematic application of the usual BFT method in principle by considering the solutions 
like as Eq.(19) [7, 13-15] although the inspection solutions may not impossible in some simple 
cases. However, in many cases like as in our non-Abelian model this calculation needs some 
routine and tedious procedure depending on the complexity of the quantities. 

In this respect, we consider our recently proposed approach using the novel property [7, 16, 
17, 20, 21] 

K[J^;^ = K[^]- (38) 

for the arbitrary function or functional K defined on the original phase space variables unless 
K has the time derivatives. Then the following relation 

{K[P],e^} = (39) 

is automatically satisfied for any function K not having the time derivatives because JF and 
their spatial derivatives already commute with 6^^ at equal times by definition. However, we 
note that this property is not simple when time derivatives exits because this problem depends 
on the definition of time derivative. This problem will be treated in Sec. IV. and we will show 
that the equality of Eq. (38) is still satisfied even with time derivatives. 

On the other hand, since the solution K of Eq. (39) is unique up to the terms proportional 
to the first class constraints 6J^, [16, 17] 

<[^]e^, (40) 

K[J-'] can be identified with K[J-'; $"] modulus the term (40). 
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Using this elegant property we can directly obtain, after some calculation, the desired first 
class Hamiltonian He corresponding to the canonical Hamiltonian H,. of Eq. (4) as follows 



ha:f-<^' 



H^m 



dPx 



iO a I ^1 a 



2«r + 4(^1)' + 2^1^"" + -'^^ 



/C 



a\2 



m 



9 



+^ M* " _ -{W'^d^wy - K + — $1 " + /c")e^ 



1 

m 



9 • 9 

/Ttl f ^ 2 ITfl 

d'x -—[{d,w-'d,wr} — Atiw-%wr 

^ lit 



where we have used the abbreviations 



/C" 



E 

n=l 



n 



1 fgV 



(n + 1)! m"^ \m 



[<|.2x($2x(...($2xf])...))];.foid 



such that Aq can be expressed as 



A'^ = A'^ + — $^ " + IC". 



m 



(41) 



(42) 



(43) 



Then, according to Eq.(3) and the property (38), the first class Hamiltonian for the total 
Hamiltonian Ht becomes 



Ht[A^^, n,; $*] = H,[A>^, n,; $*] - / d^x{d,ATQ1- 



(44) 



Note that the difference of Ht and H^ is physically unimportant since the difference is nothing 
but the ambiguity of (40) which being inherent in the definition of the BFT physical quantities 
K [16]. On the other hand, our adopted method for obtaining He using (38) as well as the 
compact forms (31)-(33) is crucial for the relatively simple and compact result which has not 
been obtained so far [8, 15]. 

Furthermore, it is important to note that all our constraints have already this property, i.e., 
e^[J^; $] = e^[^] like as follows 



0^ 



TTn 



0^-m$ 



{wnw 



fabc /\b,i;^c 



d'^t + gn^A^'^fc^ + m'A 



2 Aafi 



m 



(45) 
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In this way, the second class constraints system Oq[^] ~ is converted into the first class 
constraints one O^fJ^; $] = with the boundary conditions 0J^|$=o = 0^- 

On the other hand, in the Dirac formalism [1] one can generally make the second class 
constraint system 0a[jF] ^ into the first class constraint one 0a [J?^] = only by deforming 
the phase space {A^^, vr^) without introducing any new fields. Hence, it seems that these two 
formalisms are drastically different ones. However, remarkably the Dirac formalism can be 
easily read off from the usual BFT-formalism [7] by noting that the Poisson bracket in the 
extended phase space with $ — * limit becomes 

{A,B}\^=o = R5}-{A0a}A°"'{0a,,i?} 

= {A,B}n (46) 

where A"" = —X^'^ujpp/X^ °' is the inverse of A^a' in Eq. (6). About this remarkable relation, 
we note that this is essentially due to the Abelian conversion method of the original second 
class constraint. In this case the Poisson brackets between the constraints and the other things 
in the extended phase space are already strongly zero 

{0a,I} = 0, 
{0a, ©4 = 0, (47) 

which resembles the property of the Dirac bracket in the non-extended phase space 

{Qa,A}D = 0, 
{0a,04D = 0, (48) 

such that 

{0a,I}|$=O = {0a, A}* = 0, 
{0a,04l*=O = {0a, 0/3}* = (49) 

are satisfied for some bracket in the non-extended phase space { , }*. However, due to the 
uniqueness of the Dirac bracket [22] it is natural to expect the previous result (46) is satisfied, 
i.e., 

{ , r = { , }d (50) 

without explicit manipulation. Moreover we add that, due to similar reason, some non- Abelian 
generalization of the Abelian conversion as 

{0a, 1} = aa/S^" + aa;37$^$^ + ■ ■ ■ , 
{0a, 0/3} = /5a/37<^^ + /3a/375<^'''^' + --- (51) 

12 



also gives the same result (46) with the functions aap, ctais-y, PaP'y, etc, ■ ■ ■ of the original phase 
variables J-'. 

As an specific example, let us consider the Poisson brackets between the phase space vari- 
ables in Eq.(26). If we calculate these brackets between the BFT physical variables, after some 
manipulation we could obtain the following result 

{A'^'\x), A''^{y)} = -^{6'^'di - gr'^A^'^K^ - y). 



1 

b,k/ 



{i'^'O(x), A^'O(y)} = —gr'^ie^, - m^iS)5(x - y) 



{ir;{x)X{y)} = o, 

{I«'°(x),7f5(|/)} = -l^gr'^n^{x)6{^ - y), 

{A^^'^ix),7r',{y)} = 0. (52) 

All these results can be directly calculated by using the solutions of the BFT physical variables 
in Eq.(26). However, in some cases this is not easy work, but needs routine and tedious 
calculation as in the determination of JF in Eq.(26). This is the case of the Poisson brackets 
involving Aq, i.e., {Ag, A^'-'}, {A"'°, /1^'°}, and {A"-''^ , tt'-} . In these cases the calculations are 
more simplified by expressing A"'° as (l/m^)[02 — (-D^ttj)"] and using the Poisson brackets 
involving A"'* and tt^ which can be easily calculated from the direct use of Eq.(26). 

Here, it is interesting to note firstly that all the Poisson brackets between jF's are expressed 
still by jF's which will be found to be important later in the discussion of the time derivatives 
of the BFT physical variables. Furthermore, the results of these Poisson brackets of JF have 
exactly the same form of the Dirac's brackets between the field JF but only replacing JF instead 
of JF if there are field dependence in value of the Poisson brackets. Actually these two properties 
are general ones from the property (46), which can be checked explicitly in the case as follows 

{A'^'\x),A''\y)}\^=o = ^(5«^9^ -^r'^A^'^>(x-y) = {A^''{x),A''^{y)}o, 

{A^'\x),A'''{y)}\^=o = ^,gr'%e'2-m'A^)6{^-y) = {A^'%x),A''%y)}n, 



{A'''^ix),A'^\y)}\^=, = = {A'^'^ix),A'''\y)}n, 
{K{x)X{y)}W=o = = {n:{x),7rt{y)} 



D, 



{A'^'^(x), vf5(|/)}|^=o = 5"''5;-5(x - y) = {A«'^(a:), 7r5(i/)}z,, 



{A"'"(a:),7f;(|/)}|^=o = —,gr''7c^,{x)6{^-y) = {A^'\x),7r'^{y)}n, 
{A'^'^{x)X{y)}\^=o = = {A^'^{x),n',{y)}o. (53) 
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In this case the function K in Eqs.(38) and (39) corresponds to the Dirac brackets {A, 5}/), 
and hence K corresponding to {A, B} becomes 

{i,5} = {A5}DU^i,B^B (54) 

which proving our asserted two properties. 

Now, since in the Hamiltonian formahsm the first class constraint system without the CS hke 
term [13] indicates the presence of a local symmetry, this completes the operatorial conversion 
of the original second class system with the Hamiltonian Ht and the constraints O^ into first 
class one with the Hamiltonian Ht and the constraints O^. From Eqs. (18) and (41), one 
can easily see that the original second class constraint system is converted into the effectively 
first class one if one introduces two fields, <l>^, <l>^, which are conjugated with each other in the 
extended phase space. Note that for the Proca case the origin of the second class constraint is 
due to the explicit gauge symmetry breaking term in the action (1). 



IV. Corresponding first class Lagrangian: classical analysis 

Now, let us consider the Lagrangian (first class Lagrangian) corresponding to the first class 
Hamiltonian Ht (or Ht). It is conjectured that our first class Lagrangian must be related 
with the generalized Stiickelberg Lagrangian [19] which is gauge invariant by construction by 
noting the relation of the first class constraints and the gauge invariance (generated by the first 
class constraint) according to the Dirac conjecture [1]. To confirm this conjecture which being 
highly nontrivial in our case due to the infinitely many terms involved with the Hamiltonian 
and Lagrangian is the object of the section. 

To do this, we essentially must perform the inverse Legendre transformation by using the 
Hamilton's equations of motion as will be done in this section. But, before this traditional 
approach, we first consider the approach which showing the Lagrangian form more directly by 
considering the BFT Lagrangian formulation newly. 

A. BFT Lagrangian formulation 

In the previous sections, we have only concerned about the Hamiltonian formulation of 
the BFT method by considering the first class Hamiltonian H^ or Ht which can be easily 
calculated by replacing JF by JF according to the property (38). However we note that it is not 
clear this direct procedure of obtaining the first class quantities K[J^\ $] can be valid for the 
Lagrangian since in this case the time derivatives of the BFT physical variables are involved 
which are not clear to be the BFT physical variables also contrast to the spatial derivatives 
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of the BFT physical variables. However, interesting enough the time derivatives of the BFT 
variables remain BFT physical variables such that the procedure can be also applied to the 
Lagrangian also. 

To see this, let us consider the Poisson bracket between time derivatives of the BFT physical 
variables with the effectively first class constraints as follows 

{9°^(x),0^(y)} = I^{:F{x),H,} + J £zV{x,z)Qi{z),Qp{y)^ (/5 = 1,2) 

= {{:F{x),H,},ep{y)} 

= - {{H,,e^iy)},:Fix)} - {{epiy),:F{x)},H,} 

= 0, (55) 

where we have used the explicit form, i.e., strong form of the time derivatives as 

doH^) = {^,Ht} 

= {^{x),Hc}+ f dhV{x,z)&i{z), (56) 

V{x,z) = m^),u{z)} 

in the first line and using in the second line the fact that V should be function of only the 
BFT variables since the constituents of u and JF should be all BFT physical variables and 
no explicit JF dependence in the Poisson algebra of the fundamental BFT physical variables 
as (52). Furthermore we have also used the Jacobi identity in the third line of Eq. (55). 
However, this proof can be also performed by using the strong form of the time independence 
of the constraints (18). We note that this result that the time derivatives of the BFT physical 
variables become also the BFT physical variables is universal for any theory due to essentially 
the general validity of the formula (55). 

Now, by using this powerful property we can directly obtain the first class Lagrangian 
corresponding to the first class Hamiltonian without using the Hamilton's equations of motion. 
To this end, we first note that due to the property (55), the BFT variables ttj corresponding to 
TTj = Fio = diAo - doAi - ig[Ai, Aq] is 7r(JF)|^_^^, i.e., 

vff [^; $] = Ft, [^; $] = 9,Ig - do At + gr"^A\Al = i^« [^] . (57) 

Then, by noting the forms of (31) and (32) which looks like tiled field as gauge transformed 
one with gauge transformation matrix W, this form (57) is consistent with those forms (31) 
and (32) only when Aq is also expressed like as the gauge transformed one as follows 

I«[jr. $"] = _1 (^d^w ■ W-^y + (WAoW-^y . (58) 
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Of course, if the equality (58) is satisfied up to the constraint terms, this is also the case for 
(58). Moreover, we note that by comparing the form of (26), (58) can be explicitly checked 
by considering the time derivative of W through the Hamilton's equation of motion. Hence 
together with the forms of (31)-(33), we obtain the covariant expressions 






9 ^ 



WA^W 



-1 



(59) 
(60) 



for the four gauge fields and the field strength tensor, respectively. In this form, the gauge 
invariance of the BFT physical variables are manifest. 

Now, using this interesting result, we can consider the first class Lagrangian which is strongly 
involutive by just replacing JF by JF as follows 



C 



m 



Tr 



1 



--iWF^4J^]W"')iWF^''[J^]W-^] 



+m^ — {d^W)W-^ + WA^W'^ 
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-{d^'W)W-^ + WA^'W-^ 



Tr 



-^F^,[^]F^^[^] + m' (a, - -W-\d^W)\ (a^ - -W-^d^'W) 



, (61) 
i5«^ for 



where we have used the properties (59) and (60) and the normalization Tr{T"'T^) 

group generator T"". This provides the desired Lagrangian corresponding to the first class 

Hamiltonian He or Ht- 

On the other hand, the first class Lagrangian C can be also obtained by the inverse Legendre 
transformation of H(. or Ht as follows 



£ 



^fA - Hr 






Hn 



with the Hamiltonian density Tic and T-Ct 



(62) 



(63) 



corresponding to the Hamiltonian He and Ht of Eqs. (3) and (4), respectively. Here we have 
used (60) and the property 



d^A''^^ = 0, 



which can be obtained by replacing JF by JF in the equation 

doA''^\x) = {A"'°(x),iJT} = -d.A'^'^x). 



(64) 



(65) 
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Here, it is important to note that the inverse Legendre transformation is involved only the 
BFT variables because the first class Hamiltonian H-p or He are expressed solely by the BFT 
variables. Furthermore, we note that the first class Lagrangian C has been reduced to the 
so-called the generalized Stiickelberg Lagrangian [19] 



*~'Stuckelberg ^ ' 



--Ff.^Ff"' + m^ Ia^ - -W-\d^W)\ M^ - -W-^d^'W) 



(66) 



in the last line of Eq.(62) after in the trace operation with the identification of the field $^ 
with the Stiickelberg's scalar such that the Stiickelberg's formulation can be understood as an 
Lagrangian formulation of the BFT method. However, it seems appropriate to comment that, 
contrast to Abelian model, the simple replacement of 

A,^A^- -W-'d,W (67) 

9 

in the Lagrangian C[J-'] of Eq. (63) inspired by the mass term of (66) does not reproduce 
the generalized Stiickelberg Lagrangian since in that case the kinetic term Tr^{F^yF^^) is not 
invariant under this replacement (67) but produces several additional terms. In this sense our 
understanding the Stiickelberg Lagrangian as the original Lagrangian form in the BFT variables 
space is more simple and hence can be considered as the more fundamental formulation. 

B. Lagrangian from the Hamilton's equation of motion 

We now derive the same first class Lagrangian C from the traditional approach which is using 
the Hamilton's equations of motion [21]. However, as we shall see, some slight modification in 
the BFT Hamiltonian Ht (or H,^ should be needed in order to realize the result of (58). 

In order to understand the idea of this modification, we first consider the relatively simple 
case of Abelian theory without modification which is described by the first class Hamiltonian 



Ht = / d-^x 



1^2 , li5 2 , 1 2' r„2 -.2 



^n,' + -F,, + -m\A^ + A^ ) - A^Q^ - d^A^^i 
Then, the Hamilton's equation of motion for Ai is resulted to be 



(6J 



Ai = {Ai, Ht} = -n^ + ddo - ^-9^02 (69) 



m"^ 



such that the equality corresponding to the Abelian version of (57) is satisfied only weakly, i.e., 

n, = d,Ao - doA, - -^8,02. (70) 
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Here, the relation corresponding to (58) is not satisfied strongly but only up to non-tilled Gauss 
law constraint B2 



io = ^ + -<9o$' - ^02. (71) 



Now, in order to realize (58) in a strong form we must consider the modified Hamiltonian 

H*=Ht+ f d^x—Q2 (72) 

J m 

which has been originally introduced in other context by us recently [16]. With this Hamiltonian 
as the true time translation generator, we obtain 

A, = {Ai,H*T} = -^, + d^Ao (73) 

such that we obtain the strong equality 

TT, = 9,io - 5oA (74) 

and hence the strong equality for the gauge field Aq 

Ao = Ao + -90$' (75) 

m 

or by using the the Abelian result of Ai of (31) the the four gauge fields are found to be 

A^ = A, + -d^^\ (76) 

m 

Now, with this strong equality we can easily obtain the Stiickelberg action by considering the 
inverse Legendre transformation (62). Interestingly enough, we note that the time evolution of 
the constraints generated by H^ is isomorphic to the original ones as follows 

e2 = {e2,#*} = o. (77) 

However, we note that this modified Hamiltonian does not change the physically observable 
because the modification term affect the original Hamilton's equations of motion only by the 
terms proportional to the first class constraints also as follows 

{^(x),^^} = {^{x),Ht} + j (fy{H^),^{y)}Q2{y). (78) 

This is the idea of the Hamiltonian modification. 



r, = d,Ai - d,Ai + gr'''A\Ai - ^{D,e2r. (80) 



Now, let us return to the non-Abelian case which is described by the Hamiltonian (44) before 
modification. Similar to the Abelian case, the Hamilton's equations of motion are resulted to 
be 

if = {if, Ht} = -nt + AiS - -^iW2r (79) 

such that the equality corresponding to (57) is satisfied only weakly, i.e., 

1 
m 

Hence, in order to realize (58) in a strong form we must consider the modified Hamiltonian 

= Ht + j (fx{Al-A''^)Ql (81) 

which produces the strong equality (57) and hence our desired form (58) or (59). The Abelian 
case (72) is easily reduced form this general non-Abelian formula. Then, in this case we can 
also easily obtain the generalized Stiickelberg Lagrangian as in the previous subsection A. 
Furthermore, with this modification, we can also have the same evolution of the constraint as 
the original ones as 

§i = {ef,^*} = e^, 

e2 = {e^,^^} = o. (82) 

V. Conclusion 

In conclusion, we have applied the newly improved BFT Hamiltonian and Lagrangian for- 
malism, which provide the more simple and transparent insight to the usual BFT method, to 
the non-Abelian Proca model as an non-trivial application of our new method and completely 
solved the model by the ingenious choice of the antisymmetric X„^, u'^^ as we have suggested 
recently [16]. 

Firstly, by applying the usual BFT formalism [7] to the non-Abelian Proca model, we have 
converted the weakly second class constraint system into a strongly first class one by introducing 
new auxiliary fields. As a result we have found that the effectively first class constraints needs 
the infinite terms involving the auxiliary fields which being regular power series form in our 
ingenious choice. 
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Moreover, according to our new method, the first class Hamiltonian, which also needs 
infinite correction terms, is obtained simply by replacing the original variables in the original 
Hamiltonian with the BFT physical variables. The effectively first class constraints can be also 
understood similarly. Furthermore, we have shown that the infinite terms in the BFT physical 
variables and hence first class constraints and the first class Hamiltonian can be expressed by 
the exponential involved form. On the other hand, we have also shown that in our model the 
Poisson brackets of the BFT physical variables fields in the extended phase space are the same 
as the Dirac brackets of the phase space variables in the original second class constraint system 
only replacing the original phase variables by the BFT physical variables. We note that this 
result is conformity with the general result of BFT [7] and hence support the consistency of 
our result. 

We have also newly developed the BFT Lagrangian formulation, which has not been unclear 
so far [16, 17], to obtain the first class constraints which complements and provides much more 
transparent insight to the Hamiltonian formulation. Based on this BFT Lagrangian formulation 
we directly obtain the Stiickelberg Lagrangian as the classical Lagrangian corresponding to the 
first class Hamiltonian which is most non-trivial conjecture related to Dirac's conjecture [1] in 
our infinitely many terms involved model. This is also confirmed by considering the Hamilton's 
equations of motion with some modification term in Hamiltonian proportional to the first class 
constraints G2. This result disproves the recent argument on the in-equivalence of Stiickelberg 
formalism and BFT formalism. 

As for the other applications of our new method, we have found that our method is also 
powerful in the system with both the second and first class constraints like as the non-Abelian 
Chern-Simons gauge theory [24, 25]. Moreover, the investigation on the corresponding dual 
theory [21, 26] which would be the non-Abelian Kalb-Ramond field interacting with the Yang- 
Mills field would be interesting. Furthermore, we note that recently our improved BFT method 
is found to be also powerful in the path integral equivalence of the Maxwell-Chern-Simon theory 
with the Abelian self-dual model [27] and the BFT formulation of the non-Abelian self-dual 
model [28]. 

Finally, we would like to comment on the interesting work recently done by Banerjee and 
Barcelos-Neto [29]. After finishing our work, we have found that their work is similar to our 
work. However, we have found that their work is still incomplete compared to ours in that firstly 
they had not found our compact forms (31) and (32) and hence could not arrange explicitly the 
complicate first class Hamiltonian (41). Furthermore, they did not completely determine A^ 
contrast to our complete determination in Eq. (26). Nevertheless, it is interesting to note that 
they have succeed in proving the equivalence of BFT Hamiltonian formalism and the generalized 
Stiickelberg's Lagrangian formalism by identifying the constraints of both formalism order by 
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order. This can be considered as another independent proof of our results of section IV, where 
we have prove the result more compactly. Moreover, we have directly read off the Dirac bracket 
from the Poisson brackets of the BFT physical variables. 
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